Is the droplet theory for the Ising spin glass inconsistent with replica field theory? 
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Symmetry arguments are used to derive a set of exact identities between irreducible vertex func- 
tions for the replica symmetric field theory of the Ising spin glass in zero magnetic field. Their range 
of applicability spans from mean field to short ranged systems in physical dimensions. The replica 
symmetric theory is unstable for d > 8, just like in mean field theory. For 6 < d < 8 and d < 6 
the resummation of an infinite number of terms is necessary to settle the problem. When d < 8, 
these Ward-like identities must be used to distinguish an Almeida- Thouless line from the replica 
symmetric droplet phase. 
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Field theory has proved to be an extremely useful 
tool in studying critical transitions in ordinary systems, 
mostly by providing standard methods like the loop ex- 
pansion (above the upper critical dimension) and the 
renormalization group Its adaptation to the spin 
glass problem came just after the introduction of the 
Edwards-Anderson model and application of the replica 
trick , resulting in a kind of replica field theory (see Q 
and references therein). The first solution of the mean 
field theory of the Ising spin glass provided a simple tran- 
sition in zero external magnetic field from the paramag- 
net to a replica symmetric (RS) spin glass state |4( which, 
however, was later proved to be unstable This insta- 
bility then persisted perturbatively down to the upper 
critical dimension six and even below it 0. Repair 
came soon, at least on the mean field level, by the fa- 
mous replica symmetry breaking (RSB) scheme of Parisi 
(for details, see Q): the spin glass transition is now to a 
lower symmetry phase which is marginally stable for all 
T < T C) i.e. it is a massless phase, and has a special — 
ultrametrical — hierarchy. The RSB transition has also 
the peculiarity of extending to nonzero magnetic fields 
along the Almeida-Thouless (AT) transition line 0. 

From this point on, the spin glass community has be- 
come highly divided about the type of the transition and 
the structure of the spin glass state of short ranged, finite- 
dimensional models. Supporters of the RSB scenario fol- 
lowed the classical route trying to build a field theory on 
the basis of the — highly nontrivial — Parisi solution [(j . 



As it turned out from the physical interpretation of RSB 
, the Parisi theory corresponds to a complicated ergod- 
icity breaking with a Gibbs state decomposed into many 
pure states. On the other hand, the so called droplet the- 
ory [HQ, which was developed mostly on the original lat- 
tice system, has a simpler phase structure with only two 
pure states which are related by the spin inversion sym- 
metry (just like in a ferromagnet), and it predicts that 
a magnetic field destroys the transition. In the droplet 
picture the mode called replicon (R) remains massless in 
the whole spin glass phase |8( providing the only common 
feature both theories share. 

The simple phase structure of the droplet theory im- 
plies an RS phase with a nonzero order parameter, 
the corresponding replica field theory has a Lagrangean 
which is invariant for any permutation of the n replicas 
(up to cubic order, it was displayed in Ref. 01)- Since 
the coupling constants of such a field theory are cho- 
sen by symmetry, or — alternatively — thought to be 
the outcome of summing out short ranged fluctuations 
down to the momentum cutoff A, they have little mem- 
ory of the original control parameters like temperature 
and magnetic field. As a nonzero magnetic field does not 
change the symmetry, it is difficult to distinguish between 
a zero-field RS spin glass phase and an AT line (massless 
and replica symmetric too). In this letter, we use exact 
symmetry arguments to find out the thermal route of the 
spin glass in zero magnetic field when crossing T c from 
the paramagnet having the Lagrangean 
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where the classical fields (f) a @ are symmetrical in the replica indices a, f3 — l,...,n with zero diagonal, and 
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momentum conservation is understood in the primed 
sums. The number of spins N goes to infinity in the 
thermodynamic limit, while n — > ensures the spin 
glass limit. (A barred notation for the bare coupling 
constants is used to distinguish them from the corre- 
sponding exact vertex functions.) In addition to the 
permutational symmetry of the presumed low temper- 
ature phase, the invariants in Q have the extra at- 
tribute that a replica occurs always an even number of 
times; a Hubbard-Stratonovich-like derivation 
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makes this point evident. In mathematical form, C is 
invariant under the transformation = (— 1) Q+/3 4> a P . 
Following that transformation by the special permuta- 
tion of grouping odd and even replicas separately — i.e. 
(1, 2, 3, . . . , n) i — ► (1, 3, 5, . . . , 2, 4, 6, . . . ) — may lead us 
to figure out that we can define a class of transformations 
leaving C of invariant as follows: Let us divide the 
n replicas into two groups consisting of p and n — p el- 
ements (p being a free parameter). For the transformed 
field 6' we have: 
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for a and j3 in the same group, 
for a and (3 in different groups. 
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Therefore the paramagnetic phase has a higher sym- 
metry than even the simplest, generic replica symmet- 
ric, spin glass phase, and the presumed paramagnet to 
droplet spin glass transition breaks that higher symme- 
try in replica field theory. 

Proceeding further, we can follow the steps of Ref. 
to conclude exact identities between the irreducible ver- 
tices of the low temperature RS phase. Including a source 
term — J2 a p ^<*/3^p=o mto GJ> a Legendre-transformed 
free energy -F(q) can be derived, and it is invariant un- 
der the obviously orthogonal transformation O of Eq. (0 
1 . The derivatives of F provide the zero-momentum one- 
particle irreducible vertices Q; their definitions are the 
following: 
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and similar formulae for W a p n s,ny, U a p n s,iiv,pu: 
From -F(q') = F(q) and the orthogonality of O follows 
(using tensorial notation): 



H' = OH. 



M' = OMO 1 



(3) 



and analogous relations for higher order terms. As p is a 
free parameter, O is in fact a continuous symmetry trans- 
formation, and assuming that q is replica symmetric, it is 



1 The n(n — l)/2 q a p can be arranged into the vector q, and 
the transformation in Eq. J5J may be written then as q' = Oq 
with the diagonal tranformation matrix O having the properties 
O a f3 ia /3 = 1 or —1, depending on a and j3 belonging to the same 
group or not. 



easy to derive \J (q' — q) 2 = 2qy / p(n — p), providing an 
infinitesimal transformation for p(n — p) small. We can, 
therefore, expand the left hand sides of Eq. J2J) around 
q, and equating the coefficients of identical powers of 
p and n ~ p on both sides. Assuming stationarity, i.e. 
H = 0, and remembering that we are in a RS state with 
a nonzero order parameter q, several identities can be de- 
rived in this way between the vertices mi, m.2, wi, 
. . . , w%; ui, ... etc. 2 . These are, however, — unlike tra- 
ditional Ward-identities — power series of q with higher 
and higher order vertices, the most prominent ones are 
displayed here: 
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following from the first and 
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from the second equation in Eq. © . (Ellipsis dots in the 
above expressions are to substitute higher order terms in 
n or q.) 

The value of the identities in Eqs. Q-Q, and the 
others not displayed here, rests on their generality - 
their derivation used only symmetry arguments — ; as 
a result, they must be valid for the mean field as well 
as for low dimensional systems. It is tempting to solve 
these equations iteratively, i.e. assuming all the vertices 
are analytical in q. The most important result we can 
get in this way is the famous instability of the replicon 
mass Tr for n — * and u 2 > at criticality: 



T R = 2mi = --u 2 q 



0(q 3 )- 
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Moreover, all the vertices incompatible with the symme- 
try of the paramagnetic phase are expressible in terms of 



2 This identities are much simpler when displayed in terms of the 
set of vertices with the lower case notation (their bare counter- 
parts are the coupling constants in front of the invariants with 
the unrestricted sums in the Lagrangean). The linear relation- 
ship between these two sets of parameters were derived in 111! 
for the m's and «5's; Eqs. (20)-(24) and Table 1. As a prop- 
erty of the generic RS symmetry, it is not restricted to the bare 
couplings, but is valid for the exact vertices too. The formulae 
for the 23 quartic couplings are more complicated, and will be 
detailed in a longer publication. 
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those present for T > T c too. In leading order in q wc 
have m2 = —w\q, W2 = U2<Z, = u^q, W4 = u^q and 
W5 = Uiq, while all the others (m.3, Wq, W7, Wg) are of 
order q 2 . All these results can be verified for mean field 
theory (for a generic n) using the explicit formulae of 
Ref. [l0( and their extensions to the quartic order, and 
exploiting the fact that bare and exact parameters are 
identical for a zero loop calculation. 

We cannot follow this procedure in any finite dimen- 
sion d, as the exact vertices are no longer analytical as 
a function of q. Nevertheless, Eqs. can now be 

used perturbatively, for d > 6, to compute the bare pa- 
rameters as a function of q, and they will have, beside 
the analytical part, terms with noninteger, d-dependent 
powers. Instead of Eq. JSJ, we now have for the deviation 
of the bare replicon mass from its critical value: 
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where e = 6 — d. Cd — unlike the coefficient of the 
quadratic term — can not be computed by a simple trun- 
cation of (UJ, since a contribution from an arbitrary k- 
point vertex (multiplied by q k ~ 2 ) must be included 3 . 

At that point, two important remarks are appropri- 
ate. Firstly, there is some ambiguity in assigning the 
bare coupling constants to a given physical state below 
T c fixed by q: an offset of the zero-momentum fields, 
0p£ o — ► 0p£ o — $, leaves all the irreducible vertices 
unaltered while the bare couplings changing. We can get 
rid of "tadpole" insertions by the choice $ = q rendering 
the one-point function zero |3j. We use that case for a 
unique definition of the bare parameter space throughout 
this letter. Secondly, we must emphasize that perturba- 
tion theory remains valid for d > 6 in the loop-expansion 
sense, since Cd as well as the coefficient of the q 2 term 
in (JSJ, and in fact any quantity, can be computed, at 
least in principle, in terms of w 2 , ui,. . . ,Ui, etc. Never- 
theless, unlike the analytical contributions which can be 
computed by the truncation method just as in the mean 
field case, to get the d-dependent powers at a given order 
of the loop-expansion, we must resum an infinite number 
of terms (all the "bubbles" for Cd at 1-loop order). (We 
are not able to do this at the moment, though it may not 
be a completely hopeless task.) 

An evaluation of Eq. © at one-loop level provides us 
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with C' d again comprising an infinite number of terms, 
and there is a good chance that the linear term is exactly 
—wiq, i.e. it originates from the offset of the fields alone. 



After elucidating the g-dependence of the bare param- 
eters, expressions for the exact masses can be straight- 
forwardly computed at the one-loop level: 
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for n — and (keeping a generic n now) 



Ta — Tr = (n — 2)to 2 = — (n — 2)wiq 1 
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The one-loop fc-leg "bubble" gives a (ioi<y) 3 fc 2 contribution. 



(12) 

The dots in (|ll(l refer to the infinite terms coming from 
the 5-leg, 6-leg, etc. vertices. (Similar terms in Ijl2|l are 
not displayed, as they are definitely subleading when d > 
6.) While the anomalous mode (A) behaves regularly for 
any d > 6 and is massive, there is a competition between 
the two terms in Tr: as long as d > 8 the ci-dependent 
power in (jl 1|> is subleading, and the truncation method 
to calculate leading terms works as well as for mean field 
theory. We thus conclude that the RS phase is unstable, 
and RSB characterizes the spin glass phase for d > 8. 

Stability depends on the infinite sum in Eq. (|11|) , which 
is not available at the moment, when 6 < d < 8. We are 
taking the opportunity to comment on earlier works now. 
In Ref. 0, the second term for Tr without the infinite 
sum represented by the dots was obtained, and — as it 
is manifestly negative it was inferred that instability 
thus persisted below 8 dimensions. The authors, how- 
ever, found a sophisticated way to construct a different 
RS solution which is marginally stable, i.e. Tr = 0. (A 
modified version of that solution was proposed in a recent 
paper as a candidate model for the droplet theory (l^.l 
The first result was based on the traditional way to build 
up a symmetry-broken theory from the symmetrical La- 
grangean (Q: the bare parameters are continued analyt- 
ically below T c , and the offset of the fields, as explained 
above, then breaks the symmetry. This procedure is ob- 
viously correct for the leading terms, e.g. the 0(q) term 
of i|10|) is reproduced, and it gives definite predictions for 
the coefficients Cd, C' d , ...etc. These, however, should 
be justified by showing that they satisfy the exact identi- 
ties above. As a matter of fact, the second procedure of 
Ref. resulting in the massless replicon mode is also cor- 
rect in leading order; not surprisingly Cd is tuned to shift 
to zero. The most remarkable observation concerns 
the inclusion of an external magnetic field H cxt which 
scales as H 2 xt ~ q 2 ~^. It contributes also to Cd, and 
may render the replicon mode massless 13]. The only 
way to distinguish an AT line from a dropletlike phase is 
by means of the exact identities which exploit the extra 
symmetry of the zero-field case. 

Finally we turn now to the d < 6 case, where the loop- 
expansion breaks down, and the perturbative renormal- 
ization group elaborated in 0] takes over its role. What 
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follows from now on is highly based on the details of 
Ref . ^3 . To study the crossover region around the zero- 
field fixed point w{ 2 = — ej (n— 2), it is useful to introduce 
Wegner's nonlinear scaling fields 0] defined by the ex- 
act renormalization equations g\ — Xigi where Aj's are 
the scaling exponents. The bare parameters can be ex- 
pressed in terms of the g^s, as displayed here for the three 
masses (omitting nonlinear terms and irelevant fields): 

?1 — 2 _ 2 

mi = w 1 + (gi + 2g 2 + g 3 j + ■■■ , 

m 2 = -(52+53) + • • ■ , 
1 

m 3 = ^93 + ■ ■ ■ ■ 

Assuming the form gi = Ci(w\q) Zi , we can conclude from 
(JSJ) that C3 = O(e), hence hindering us to compute the 
correction term of 23 = 1 + . . . at that order. The other 
two exponents, however, can be derived from the loga- 
rithms of @ and J5j) providing Z\ = 1 — e/2 + 0(e 2 ) and 
22 = 1 + 0(e 2 ). These results are in accord with the re- 
lations A, = (2 — e/2 + ri/2) z% — which must follow from 
the flow of the order parameter q = ((3/v) q, with (3/v = 
(2 — e/2 + 77/2) — when comparingwith the independent 
calculations of the A.j's and 7/ in 14]. Although the lead- 
ing terms for the C,'s are trivial [C% = (n — 2) + O(e), 
C2 = 1 + 0(e) and C3 = 0(e)], the 0(e) corrections can 
be obtained again, just like for 6 < d < 8, including all 
the terms of Eqs. To see this, let us consider the 

scaling form of a generic fc-point vertex 

r « = | ffl |£[*-fc(2-m)] rM(x, y ), 

where we have used the common notation T for the 
vertices (instead of the m's, w's, it's, etc). The two 

scaling variables are x — (te/lffil* 1 = ffa/lsir 1 an d 

V = ,93/l5i|~ — 93/\9i\~ 1 whereas the f( fc )'s are scal- 
ing functions specific to the given vertex from the fc-point 
family. As we have for the leading terms in e: ~ ?D* fc 
(k > 3 and forgetting now about the only exception u>i), 
a typical term of, say, Eq. (QJ : 

i.e. all terms are the same e order. 

To test stability below six dimensions, Tji(x.y) must 
be computed in e-expansion first, and then substituting 
x = C2/|Ci|*i, y — Cz/\Ci\^ provides us the replicon 
mass for T < T c in zero magnetic field. From x = 1/(2 — 
n) + 0(e) and y = O(e) we can conclude: 

~ TL 

r R = (-1 + 2x + y) + ■ ■ ■ = + 0(e), 

2 — n 

In the spin glass limit n — * 0, stability depends on the 
correction term which contains two contributions: the 
first one comes from the correction of the scaling func- 
tion Tr, — which has been computed and gives a negative 



(unstable) result — , and the second one from those of x 
and y. These are, however, — as argued above — not 
available at the moment, as they result from a resum- 
mation of an infinite number of terms in the identity of 
Eq. 0}. 

To conclude, the importance to find the correct tra- 
jectory in the bare parameter space in absence of an ex- 
ternal magnetic field (what is called here the thermal 
route) is emphasized when the stability of the low tem- 
perature RS phase — the droplet phase — is tested. We 
argued that for 6 < d < 8 and d ^ 6 an infinite num- 
ber of one-loop graphs should be resummed to settle the 
problem. Without using the Ward-like identities derived 
in this letter, the droplet phase cannot be distinguished 
from an AT line for d < 8. Our conclusions are in con- 
flict with those of a recent paper by M. A. Moore [m ]. 
Although both of us realized the necessity to resum an 
infinite number of terms, in Ref. |l6| self-energy graphs 
with different loops are proposed to be summed to get 
rid of infrared divergences caused by the replicon mode. 
In this letter, we argue that allowing for all the one-loop 
graphs with different number of legs is a must to ensure 
the extra symmetry imposed by the lack of an external 
magnetic field. Old beliefs that replica field theory is 
inconsistent with the droplet picture were based on the 
instability emerging at one-loop order if Eq. Ijlljl without 
the infinite terms is used, which we believe is not correct. 
Therefore the conclusion that the RS phase is unstable 
for 6 < d < 8 is thus premature, just as that it is stable 
for d < 6. Eq. Hll|) loses its relevance for d < 6, and 
the study of the crossover region around the nontrivial 
fixed point is inevitable. It must be stressed that our 
identities are valid even in 3 dimensions, though pertur- 
bative methods are not available then to take advantage 
of them. 

Discussions with C. De Dominicis are highly appre- 
ciated. I am grateful to M. A. Moore for sending me 
comments to the e-print version of that letter. 
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